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摘要 :目前 对 火星 CH4 气 体 的 探测 是 探索 火星 生命 的 一 项 重要 手段 ， 圈 定 出 火星 表面 CHs 源 区 的 


位 置 可 为 将 来 火星 生命 的 探索 选取 合适 的 目标 点 。 本 文 在 对 火星 CH4 气 体 共振 散射 进行 探测 的 基 


础 上 ， 通 过 数值 模拟 的 方法 对 火星 CH4 气 体 的 空间 分 布 进行 了 反 演 。 反 演 结 果 能 再 现 模型 的 密度 
分 布 ， 辨 认 出 CH4 气 体 密度 分 布 较为 稠密 的 区 域 ， 从 而 可 以 确定 出 火星 表面 CH4 源 区 的 位 置 。 
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人 类 一 直 在 试图 寻找 地 球 以 外 是 否 有 生命 的 存在 。2008 年 , 美国 “凤凰 号 ”火星 着 陆 器 在 火 


曾经 存在 过 生命 将 是 突破 性 的 重大 发 现 。 


星 表面 发 现 冰 水 ， 找 到 了 生命 存在 的 条 件 中 。 在 未 来 的 火星 探测 任务 中 ， 寻 找 火星 上 是 否 存在 或 


在 地 球 上 ，90% 的 CH 气体 是 生物 体 在 消化 养分 的 过 程 中 释放 出 来 的 。 火 星 大 气 中 已 经 确认 


存在 CH 气体 ， 但 不 能 确认 是 否 来 自生 命 活 动 趾 。 因 此 ， 火 星 上 CH 气体 的 探测 成 为 探索 火星 生 


命 的 一 项 重要 手段 。 
前 


的 探测 结果 ， 火 星 大 气 中 CH4 的 含量 约 为 十 亿 分 之 一 品 ， 但 该 结果 仅 为 全 球 性 的 分 


布 情况 。 如 果 我 们 知道 CH4 的 全 球 空间 密度 分 布 情况 ， 就 可 以 找到 CH4 最 为 稠密 的 区 域 ， 这 样 就 
可 以 圈定 出 火星 表面 CHs 源 区 的 位 置 ， 为 将 来 探索 火星 生命 选取 合适 的 目标 点 。 从 而 为 进一步 探 


索 火 星 生 命 研究 提供 重要 的 科学 依据 。 


在 深 空 探测 领域 , 利用 气体 分 子 在 受到 辐射 后 产生 的 共振 散射 特征 谱 线 进行 探测 有 着 悠久 的 


历史 外。Drossart 等 人 在 对 木星 和 土星 的 光谱 进行 研究 时 发 现 了 CH4 分 子 荧光 光谱 的 存在 。1995 


年 ， 欧 空 局 发 射 了 红外 空间 天 文 台 太空 望 眼镜 ， 该 卫星 首次 成 功 对 CH4 红 外 波段 3.3 
行 了 探测 趾 。 目 前 ， 本 项 目 组 正 试图 研究 如 何 对 火星 CHs 气 体 的 共振 散射 进行 探测 ， 
正 是 在 此 基础 上 展开 的 。 


hm 的 荧光 进 
本 文 的 工作 


由 于 火星 CH 气体 比 较 稀薄 ， 且 不 会 发 生 二 次 散射 ， 因 此 可 近似 为 光学 薄 层 ， 探 测 器 获得 的 
数据 是 对 应 探测 方向 上 所 有 CH4 共 振 散射 强度 的 线 积 分 ， 正 比 于 相应 方向 上 CH4 的 柱 密度 。 可 以 


表示 为 ， 


P(D = af n(x)dx (1) 
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其 中 ，P(D 为 探测 器 沿 [ 方 向 探测 得 到 的 数据 ，m(x) 为 空间 x 处 CH4 的 密度 ，o 为 比例 因子 。 

这 个 过 程 和 发 射 式 计算 机 断层 成 像 ECT (Emission Computed Tomography) 的 投影 过 程 一 致 ， 
当 卫 星 从 不 同 角度 对 火星 CH4 进 行 探 测 时 ， 可 以 通过 CT (Computed Tomography) 方法 反 演 出 CH4 
的 空间 分 布 。 

本 文采 用 数值 模拟 的 方法 ， 利 用 CT 方法 反 演 CH 的 空间 密度 分 布 ， 并 与 预 设 的 CH 密度 模型 
进行 比较 ， 从 而 为 观测 方案 提供 理论 依据 。 

1 CHARA! 

由 于 缺乏 详实 的 观测 数据 ， 目 前 任何 火星 大 气 模型 中 都 不 包含 CH4， 因 此 我 们 需要 预 设 出 一 
个 CH 分 布 模型 。 本 文 的 模型 是 根据 火星 大 气 分 布 规律 建立 的 ， 不 一 定 和 实际 的 CH4 空 间 分 布 完 
全 相符 。 但 事实 上 ， 在 CT 重建 领域 ， 普 遍 使 用 仿真 模型 来 研究 反 演 的 精度 。 这 些 仿真 模型 一 般 
采用 和 研究 目标 具有 相近 密度 的 材料 来 替代 , 按照 研究 目标 的 结构 进行 设计 和 加 工 而 成 , 也 就 是 
说 算法 的 有 效 性 并 不 会 随 模型 的 变化 而 发 生 本 质 的 改变 。 换 名 话说 , 如 果 将 来 卫星 采集 到 的 数据 
量 和 文中 模拟 所 使 用 的 数据 量 一 样 , 那么 反 演 结果 的 精度 也 和 文中 采用 模型 反 演 结果 的 精度 一 样 。 

我 们 采用 MCD(Mars Climate Database) 数 据 库 中 的 火星 大 气 分 布 数据 。MCD 4.3 包 含 从 火星 
表面 到 250km 以 上 的 气象 数据 ， 包 括 : URES UG. Bo. SRK CRR, KK 
汽 、 含 冰 量 ) 。 该 数据 库 的 数据 来 自 火星 大 气 GCM (大 气 环流 ) 模型 运行 的 结果 ， 与 观测 数据 
非常 吻合 ， 可 以 精确 地 给 出 不 同位 置 火星 大 气 CO,、N，,，、O、CO 体 积 比 的 廊 线 。 利 用 MCD 先 计 
算 各 经 纬度 上 大 气 的 总 质量 密度 CO，,、N，,、O、CO、0O;、 冰 、 水 汽 的 体积 比 随 高 度 的 分 布 。 其中， 
体积 比 = 某 种 成 分 的 疏 密度 /大 和 气 总 数 密 度 。 在 上 述 大 气 数 据 的 基础 上 ， 作 以 下 假设 ; 首先 假设 
CH4 非 全 球 均匀 分 布 ， 沿 着 高 度 方向 体积 比 不 变 ， 均 为 10ppb。CH4 的 体积 比 随 经 、 纬 度 为 麦克 斯 
圳 分 布 ， 其 次 ，CH4 在 大 气 中 所 占 比 例 非常 小 ， 忽 略 增加 CH 而 引起 的 大 气 总 质量 密度 的 变化 。 
在 这 样 的 假设 下 , 先 由 总 质量 密度 计算 出 大 气 的 总 数 密度 , 然后 再 按照 CH 的 体积 比 , 可 得 到 CH4 
在 不 同 经 度 、 纬 度 、 高 度 的 数 密度 。 

根据 现 有 的 观测 结果 ， 假 定 CH4 分 布 主要 集中 在 低 纬 度 的 区 域 。 因 此 模型 中 CH4 分 布 只 选取 
[50S, 50N] 的 范围 。 其 次 ， 由 于 卫星 轨道 选择 晨 昏 面 上 的 极 轨道 〈 见 3.1 节 ) ， 随 着 火星 的 自转 ， 
所 有 经 度 上 的 反 演 精度 都 是 一 样 的 ， 为 了 减 小 计算 量 ， 模 型 中 CH4 分 布 只 选取 [0E, 120E]， 并 以 
(ON, 60E) 为 中 心 。 另 外 ， 随 着 气压 降低 ，CH4 密 度 是 在 降低 的 ， 反 演 的 目的 也 是 为 了 获取 其 扩散 
过 程 ， 并 最 终 反 演 得 到 气体 源 区 ， 即 密度 分 布 最 为 稠密 的 区 域 。 为 了 降低 计算 量 ， 模 型 中 CH4 
分 布 的 高 度 极限 于 从 火星 表面 到 150km 的 范围 。 如 图 1 所 示 〈 取 自然 对 数 的 结果 )， 其 密度 的 最 大 
值 在 火星 表面 模型 中 心 附近 的 (IN,56E) 处 ， 密 度 值 为 1.9 x 1016m. 


图 1 KE CH 大 气 的 分 布 模型 

2 反 演 方法 
严格 来 讲 ， 要 做 精确 的 反 演 ， 要 求 对 经 过 待 反 演 物体 的 任意 点 的 投影 线 必须 覆盖 180。 的 角 
范围 ( 称 为 投影 数据 完备 性 条 件 ) 1， 但 通常 情况 都 不 能 满足 这 一 条 件 ， 因 此 无 法 获得 准确 的 


反 演 结果 。 但 是 ， 反 演 过 程 中 可 采取 一 些 补偿 方式 ， 通 过 某 些 先 验 知识 来 改善 反 演 的 结果 ， 提 高 
反 演 的 准确 性 。 
CT 理论 经 过 几 十 年 的 发 展 已 经 趋 于 成 熟 ， 并 成 功 应 用 于 各 个 领域 。 对 于 完备 的 投影 数据 ， 
可 以 利用 精确 地 反 演 算法 ， 如 PI 线 四、 滤波 反 投 影 算 法 〈(FBP) 巴 等 等 。 如 果 投 影 数据 有 缺失 ， 
则 通常 用 迭代 算法 来 解决 ， 如 代数 迭代 法 (ART) 钙 芍 、 向 前 投影 法 (GP) 8239 等 。2006 年 ， 
Sidky 等 人 [将 图 像 总 变 差 (TV) 极 小 化 与 迭代 算法 中 较为 广泛 应 用 的 ART 算 法 相 结合 ， 在 投影 
数据 严重 缺失 的 CT 反 演 上 取得 了 极 佳 的 效果 。 卫 星 在 观测 火星 CH4 气 体 的 过 程 中 ， 很 难 采集 到 完 
备 的 投影 数据 量 ， 因 此 本 文 将 采用 该 方法 进行 数值 模拟 。 
2.1ART 算 法 
ART 算 法 的 本 质 是 解 线 性 方程 组 。 当 探测 器 在 空间 某 一 位 置 i 对 待 重建 物体 进行 探测 时 ， 每 
一 个 探测 单元 的 探测 值 和 被 反 演 的 区 域 可 表示 为 : 
pii = hjin (2) 
其 中 ，n 为 待 反 演 区 域 体 素 密度 值 的 向 量 ，pji; 是 对 第 j 号 探测 器 的 探测 值 ，hij; 为 相应 的 投影 
和 矩阵 向 量 。 联 立 所 有 位 置 的 探测 结果 ， 可 以 得 到 ; 
P=Hn (3) 
ART 的 特点 是 先 假设 一 个 初始 值 n(%), 然后 根据 该 初始 值 求 出 第 一 次 近似 值 nz) ， 再 根据 mn 
求 出 第 二 次 近似 值 n@, 如 此 继续 ， 直 到 满足 收敛 条 件 为 止 。 在 根据 nw-? 求 %(9 时 ， 利 用 下 面 的 
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ates, 


nk = nk? =0 
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nk = nk-1 + AK 1i H HT IH, ll +0 (4) 
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2.2 TV 方法 

TV 最 早 由 Rudin 等 人 提出 来 ， 用 于 图 像 去 品 ， 后 来 逐渐 成 为 一 种 广泛 的 优化 准则 ， 用 于 提高 
CT 图 像 反 演 的 图 像 质量 。 

将 TV 约束 引入 ART 算 法 的 过 程 中 ， 也 就 是 在 投影 方程 的 约束 条 件 下 ， 求 解 图 像 总 变 差 最 小 


的 优化 问题 ， 即 


minp=np nllnllrv (5) 
其 中 ， 
lInl|lry = ys |Vn| = > ome E xyz)2 T (Ny,z = xy-1z)2 + (Ny yz aj Nyy 2-1)” 
XYZ XYZ 


通过 在 ART 算 法 过 程 中 添加 先 验 条 件 〈 在 火星 CH4 气 体 的 反 演 中 , 先 验 条 件 为 : 密度 值 非 负 ; 
若菜 一 投影 射线 上 的 值 为 0， 则 经 过 该 射线 的 所 有 体 素 的 密度 值 均 为 0) ， 引 入 TV 约束 ， 得 到 投 
影 数 据 不 完备 情况 下 更 好 的 重建 图 像 。 


3 计算 参数 的 选取 


在 卫星 轨道 和 探测 器 的 设计 方面 ， 项 目 组 采用 以 下 方案 。 
3.1 卫星 轨道 参数 
在 模拟 中 ，CH4 气 体 被 认为 随 着 火星 一 起 自转 ， 并 且 在 几 天 内 其 密度 分 布 没 有 变化 。 因 此 火 
星 同 步 轨道 卫星 无 法 观测 到 火星 表面 所 有 区 域 的 CH4。 如 果 采 用 极 轨道 ， 随 着 火星 的 自转 ， 便 可 
以 观测 到 所 有 区 域 的 CHe。 为 了 简化 计算 ， 轨 道 形状 采用 圆 轨道 。 为 了 方便 描述 卫星 所 运行 的 轨 
道 ， 定 义 一 个 火星 太阳 坐标 系 MS): 坐标 原点 为 火星 中 心 ，Xws 轴 指向 太阳 ，Yws 轴 指向 火星 
卉 人 出 ，Zws 轴 按 右手 规则 给 出 。 为 了 避免 太阳 光 直 射 探测 器 , 将 卫星 运行 轨道 设 在 晨 错 Ys-Zms 平 
面 上 。 在 模型 中 ，CH4 气 体 分 布 的 最 大 高 度 为 169km。 为 了 保证 在 轨道 的 任意 一 处 ， 探 测 器 的 每 
个 探测 单元 都 可 以 探测 到 对 应 视线 方向 上 所 有 的 CH4， 卫 星 的 轨道 高 度 设 定 为 500km， 轨 道 周 期 
1.5h. 
3.2 探测 器 参数 

模拟 中 ， 探 测 器 的 张 角 采 用 135 ， 这 样 才能 保证 在 轨道 上 任意 一 处 都 可 以 探测 到 所 有 CH4 气 


> 体 ,除了 被 火星 遮挡 的 区 域 (火星 遮挡 分 为 直接 谈 挡 和 间接 遮挡 。 直 接 谈 挡 是 ， 由 于 探测 器 探测 
© 到 的 数据 是 对 应 探测 方向 上 所 有 CH 共振 散射 强度 的 线 积分 。 那 么 ， 一 旦 探测 器 的 视线 受到 火星 
= 的 遮挡 ， 散 射 强度 线 积分 的 终点 将 是 火星 表面 靠近 探测 器 的 一 侧 ， 另 一 侧 的 CH 无 法 被 探测 到 。 
© 间接 遮挡 是 ， 在 火星 的 背面 ，CH4 不 会 与 太阳 光 产 生 共振 散射 作用 ， 因 此 探测 器 无 法 探测 到 火星 
pt 背面 的 CH4) 。 每 个 探测 单元 的 张 角 为 3 ( 即 探测 器 的 水 平分 辨 率 和 垂直 分 辩 率 均 为 50km， 这 是 
~ 项 目 组 目前 研制 的 探测 器 所 能 达到 的 分 辨 率 ), 则 探测 器 共有 45x 45 个 探测 单元 。 在 探测 过 程 中 ， 

探测 器 的 中 心 指向 火星 的 中 心 ， 采样 时 间 间 隔 为 .5min, 采集 数据 的 总 时 间 为 火星 自转 一 周 ( 约 


24h)， 因 为 在 相对 火星 不 动 的 坐标 系 下 ， 卫 星 的 轨迹 以 24h 为 周期 重复 运动 ， 如 图 2 所 示 。 


图 2 卫星 运行 轨迹 


4 数值 模拟 结果 


反 演 的 空间 范围 完全 覆盖 了 模型 中 CH4 的 分 布 区 域 ， 从 火星 表面 3379.720Kkm 到 3549.720km 
的 范围 。 反 演 中 ， 经 度 和 纬度 的 网 格 划分 均 为 每 隔 1"， 高 度 为 每 10km 一 层 ， 共 18 层 , 网 格 的 数量 
为 18 x 180 x 360。 反 演 的 结果 如 图 3 所 示 。 可 以 看 出 , 反 演 出 来 的 CH4 主 要 集中 在 纬度 [50S, 50N]、 
经 度 [0E,120E] 的 范围 内 ， 与 原 模型 基本 相似 。 


图 3 火星 CHa 大气 的 反 演 结果 

为 了 进一步 对 比 反 演 的 结果 ， 选 取 纬 度 为 0 和 经 度 为 60" 这 两 个 截面 来 分 析 。 分 别 如 图 4 上 、 
下 所 示 。 其 中 ， 左 列 为 模型 ， 右 列 为 反 演 结果 。 从 图 可 进一步 清楚 的 看 出 ，CH4 基 本 上 都 在 纬度 
[50S,50N]， 经 度 [0E, 120E] 的 范围 内 。 另 外 ， 与 模型 一 样 ， 在 相同 的 经 纬度 上 密度 随 着 高 度 的 增 
加 逐渐 降低 。 但 反 演 结果 的 密度 偏 大 , 特别 是 在 离 火 星 表 面 较 远 的 地 方 。 这 主要 是 投影 数据 偏 少 
造成 的 。 从 图 2 卫星 运行 的 轨迹 可 以 看 出 ， 卫 星 轨 道 在 经 度 方 面 有 较 大 的 间隔 ， 在 并 道 面 高 达 
11.25*， 这 导致 了 观测 数据 有 较 大 的 和 欠缺。 另外， 火星 的 遮挡 也 是 导致 投影 数据 缺失 的 另 一 个 重 
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图 4 CH 密度 分 布 在 纬度 为 0 CERD 和 经 度 为 60° 的 截面 (下 图 ) 

( 左 图 为 模型 ， 右 图 为 反 演 结果 ) 
如 果 相 机 采样 的 时 间 间 隔 可 以 缩小 , 我 们 便 可 以 得 到 更 多 的 投影 数据 , 反 演 结果 的 精度 也 会 
相应 提高 。 当 采样 的 时 间 间 隔 减 小 一 半 ， 即 从 1.5min 减 小 到 0.75min 时 ， 若 用 平均 误差 ( 即 反 演 
结果 和 原始 数据 差 值 的 点 平均 值 占 原始 数据 点 平均 值 的 比重 ) 来 衡量 反 演 的 精度 , 则 采样 间隔 为 
1.5min 的 反 演 误差 为 16.95%， 而 0.75min 采 样 间隔 的 反 演 误差 为 15.07%。 另 外 ， 如 果 考 虑 卫星 
的 进 动 ， 同 样 也 可 以 提高 反 演 的 精度 。 在 计算 中 ,我们 简单 地 假设 卫星 每 天 进 动 3"。 由 于 卫星 的 
进 动 ， 在 相对 火星 静止 的 坐标 系 中 ， 卫 星 的 运动 轨迹 不 发 生 重 合 。 考 虑 到 CH4 在 空间 上 不 断 的 扩 
散 ， 因 此 ， 采 样 的 时 间 不 宜 太 长 ， 这 里 只 采集 了 3 天 的 数据 量 ， 轨 道 倾角 从 85° 逐 渐变 化 到 94°。 
当 采 样 间隔 分 别 取 1.5min 和 0.75min 时 ， 反 演 结果 的 平均 误差 分 别 为 15.89% 和 14.34%。 表 1 即 为 
这 4 种 反 演 的 对 比 ， 从 表 中 可 看 出 采用 更 多 数据 量 ， 更 小 的 采样 间隔 ， 或 者 采用 进 动 会 使 反 演 效 
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#1 不 同 参数 反 演 的 结 


参数 1 2 3 4 
轨道 倾角 (9 90 90 85-94 85-94 
采样 数据 量 (d) 1 1 3 3 
采样 间隔 (min) 1.5 0.75 1.5 0.75 
相机 张 角 (9 135 135 135 135 
Pua AVEC) 0 0 3 3 
平均 误差 (9%) 16.95 15.07 15.89 14.34 


5 结论 和 讨论 
本 文 利用 CT 方法 对 火星 CH4 气 体 空间 分 布 的 反 演进 行 了 数值 模拟 。 数 值 模拟 结果 表明 ， 当 卫 
星 运行 于 火星 晨 昏 面 的 圆 轨道 时 , 反 演 的 密度 分 布 和 模型 基本 一 致 , 尤其 是 密度 变化 趋势 是 一 致 
的 。 根 据 反 演 的 结果 ， 可 以 明显 看 出 CH4 气 体 密度 分 布 较为 稠密 的 区 域 ， 这 可 以 满足 “圈定 CH4 
的 源 区 位 置 ， 为 将 来 探索 火星 生命 选取 合适 的 目标 点 ”这 一 科学 任务 。 
旦 是 反 演 仍然 存在 一 定 的 问题 。 首先 , 反 演 结果 存在 一 定 的 偏差 。 只 有 增加 采样 的 数据 量 ( 文 
中 考虑 了 缩小 采样 的 时 间 间 隔 ， 改 变 卫星 轨道 ) 才 可 以 提高 反 演 的 精度 。 其 次 ， 反 演 结果 的 空间 
分 辨 率 还 有 待 于 提高 。 反 演 结果 的 空间 分 辨 率 受 控 于 探测 器 的 空间 分 辨 率 , 而 本 文 是 鉴于 项 目 组 
目前 探测 器 设计 所 能 达到 分 辨 紊 的 基础 上 进行 模拟 的 (目前 , 项目 组 所 研制 的 探测 器 所 能 达到 的 
最 高 水 平分 辨 率 和 垂直 分 辨 率 均 为 50km， 如 何 提高 探测 器 的 分 辨 率 问 题 是 项 目 组 正在 解决 的 难 
题 ) ， 所 以 空间 分 辨 率 不 高 。 再 次 ， 文 中 只 考虑 了 稳 态 的 情况 ， 但 现实 中 ， 火 星 表面 的 中 性 风 会 
影响 CH4 的 空间 分 布 ， 从 而 导致 CH4 的 分 布 随时 间 变 化 。 项 目 组 将 会 在 下 一 步 研究 工作 中 进一步 
研究 这 些 问题 ， 以 期 达到 最 好 的 重建 效果 。 
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